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Abstract 



We present predictions for the inclusive cross section for Higgs boson production by 
gluon-gluon fusion in proton collisions at y/s = 8 TeV. Our calculation is accurate up 
to next-to-next-to-leading order in QCD perturbation theory and includes soft-gluon 
effects up to next-to-next-to-leading logarithmic accuracy and two-loop electroweak 
corrections. The dependence on heavy-quark masses is taken into account exactly up 
to next-to-leading order and next-to-leading logarithmic accuracy, and a treatment of 
the Higgs boson line-shape is provided according to the complex-pole scheme. 
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The data collected at the LHC in 2011 allowed the ATLAS pQ and CMS [2 J collaborations to 
impose strong constraints on the allowed mass of the Standard Model (SM) Higgs [3] boson, by 
essentially excluding it in the mass range (9(130 GeV) < tuh < 0(600 GeV). Both collaborations 
observed an excess of events around tuh ~ 125 GeV and the update of the Tevatron [4] results 
with up to 10 fb _1 integrated luminosity also points to an excess in the region 115 — 135 GeV. 

With the LHC running at 8 TeV and expecting to deliver an integrated luminosity of (9(15 ftT 1 ) 
per experiment, the discovery or exclusion of the SM Higgs boson can really be envisaged by the 
end of 2012. 

The dominant mechanism for SM Higgs boson production at hadron colliders is gluon-gluon 
fusion [5], through a heavy-quark (mainly, top-quark) loop. The QCD radiative corrections to the 
total cross section have been computed at the next-to-leading order (NLO) in Refs. [SIEIIH] and at 
the next-to-next-to-leading order (NNLO accuracy) in [HI QUI EH]- NNLO results at the exclusive 
level can be found in Refs. p2J M, M CES] . 

In this article we present state of the art predictions for this channel with explicit results 
at the LHC with centre-of-mass energy of 8 TeV. Other recent results for the Higgs production 
cross section at the LHC, obtained at different levels of theoretical accuracy, are presented in 

Refs. ma mi nans]. 

The main features of our calculation have already been described in Ref. [2UJ. Here we focus 
on the improvements with respect to the work of Ref. [20] . 

Our calculation is based on the resummation of soft-gluon contributions to all orders, as a 
way to improve state of the art fixed-order predictions with the dominant effect from higher-order 
corrections. The resummation of soft-gluon effects is achieved (see Ref. [21J for more details) by 
organizing the partonic coefficient function in the gg channel in Mellin spac^j] as 

G mn(<X8(j4i)i m \l^ WhIPf) = a s(/4) CWOsOfl), m 2 H /fi 2 R ; m 2 H /n F ) 

■ A%(a s (» 2 R ),m 2 H /fam 2 H /» 2 F ) +0(1/N) , (1) 

where as is the QCD coupling evaluated at the renormalization scale fi R and fip is the 
factorization scale. The large logarithmic corrections (that appear as a§ In" N in Mellin space) 
are exponentiated in the Sudakov radiative factor A^, which depends only on the dynamics of soft 
gluon emission from the initial state partons. The hard coefficient C gg includes iV— independent 
terms arising from both soft and hard gluon emission and depends on the details of the coupling 
to the Higgs boson and, therefore, on the masses of the heavy quarks in the loop. The coefficient 
needed to perform the calculation up to NTL can be obtained from the corresponding fixed order 
computation to N l LO accuracy. 

The explicit expression for the coefficient C gg computed in the large-m t limit up to two-loop 
accuracy can be found in Ref. [2T] . Together with the knowledge of the radiative factor A^, this 
result allowed us to perform the resummation up to next-to-next-to-leading logarithmic (NNLL) 
accuracy and to match the ensuing result with the fixed-order NNLO computation [9j [TOj [11] 
still performed in the large-m 4 limit [21]. The resummation effect has been confirmed by the 

^The Mellin transform is denned with respect to the variable z = mj^/s, s being the partonic centre-of-mass 
energy. 
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computation of the soft-gluon terms at N 3 LO [22]. The NNLL+NNLO result of Ref. [2T] has been 
the reference theoretical prediction for the gluon fusion cross section for few years. 



In Ref. [20] we have extended this calculation by including the bottom-quark loop (and its 
interference with the top contribution) exactly up to NLO [7, 8J. However, the NLO results in 
[HI EHJ EU |25] can be used to extract the exact expression of C gg in Eq. (JTJ) up to NLL: 

c 99 (as(/4),™y/4; m V/4) = 



m, 



1 + Z*™. I C , K ) + 6C(2) + i n ™ + 67 £ + ^ _ e 7s i n ^ + 0(q 2) (2) 

where 7# ~ 0.577216 is the Euler's constant and the dependence on the heavy-quark masses 
appears in the function c^(m q ). The corresponding expression can be found in Eq. (B.2) of 
Ref. [8] in terms of one-dimensional integrals, or as a fully analytic result in Eq. (3.5) of Ref. [23] 
and Eq. (27) of Ref. [24], both in terms of harmonic polylogarithms. In the limit of infinite quark 
mass, one recovers the well know result 

c^rrig) > ^ . (3) 

m q — >oo Z 

By using Eq. (J2J we can achieve NLL+NLO accuracy without relying on the large-m t approxi- 
mation, that is, we can improve the exact fixed order NLO result by including soft-gluon resum- 
mation at NLL. Since no exact results are available at NNLO accuracy, at this order we only 
include the top-quark contribution in the mt — > 00 limit by adding soft-gluon effects at NNLL, as 
in Refs. [2TJ [2U] . The ensuing result is finally corrected for two- loop electro- weak (EW) contribu- 
tions [26] as evaluated in [27], in the complete factorization scheme, in which the EW corrections 
are applied to the full QCD corrected cross section^. We point out that the inclusion of the exact 
dependence on the top- and bottom-quark masses up to NLL accuracy results in a decrease of 
the cross section ranging from about 1.5% at mu = 125 GeV, to about 6% at tuh = 800 GeV. 
The usually neglected charm-quark contribution, which we also include for the first time in our 
calculation, further decreases the cross section by about 1% for a light Higgs, being instead very 
small in the high-mass region. 

The second improvement with respect to the work of Ref. [20] regards the treatment of the 
Higgs boson width. While the Zero Width Approximation (ZWA) can be considered sufficiently 
accurate for the evaluation of the inclusive cross section for a light Higgs boson, the increase 
of the Higgs boson width at large masses requires a proper inplementation of the corresponding 
line-shape. In this work we rely on the OFFP scheme described in Ref. [2_H] as an effective 
implementation of the complex-pole scheme. The signal cross section can be written as 

/ \ f ,,,2 MT h(M) a(M) 

where mn+i^H parametrizes the complex pole of an unstable particle, with m# an input parameter 
playing the role of the on-shell mass while jjg, equivalent to the on-shell width, is computed at one 
loop accuracy in the SM in Ref. [29]. r#(M) and a(M) correspond to the conventional on-shell 
width and hadronic production cross section evaluated at the virtuality of the Higgs boson M, 
respectively. The calculation in Ref. [29] provides a realistic estimate of the complex-pole width 7^ 



■'"Results including this improvement for y/s — 7 and 14 TeV have been presented in Ref. 
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above the ZZ threshold but might introduce an artificial effect at low masses due to the fact that 
in general ~£{t(Ih)Hh 7^ 1- hi order to recover the ZWA for light Higgs we use an extrapolation 
of the value of 7^ towards the on-shell decay width T{m H ) below m# = 200 GeV [30] [§. 



m H (GeV) 


a (pb) 


scale(%) 


PDF+a s (%) 


115.0 


22.68 


+7.4 


-8.1 


+7.6 


-6.8 


116.0 


22.31 


+7.4 


-8.0 


+7.5 


-6.8 


117.0 


21.93 


+7.4 


-8.0 


+7.5 


-6.8 


118.0 


21.58 


+7.4 


-8.0 


+7.5 


-6.8 


119.0 


21.23 


+7.3 


-8.0 


+7.5 


-6.8 


120.0 


20.88 


+7.3 


-7.9 


+7.5 


-6.9 


120.5 


20.72 


+7.3 


-7.9 


+7.5 


-6.9 


121.0 


20.56 


+7.3 


-7.9 


+7.5 


-6.9 


121.5 


20.39 


+7.3 


-7.9 


+7.5 


-6.9 


122.0 


20.24 


+7.3 


-7.9 


+7.5 


-6.9 


122.5 


20.08 


+7.2 


-7.9 


+7.5 


-6.9 


123.0 


19.92 


+7.2 


-7.9 


+7.5 


-6.9 


123.5 


19.76 


+7.2 


-7.9 


+7.5 


-6.9 


124.0 


19.61 


+7.2 


-7.9 


+7.5 


-6.9 


124.5 


19.46 


+7.2 


-7.9 


+7.5 


-6.9 


125.0 


19.31 


+7.2 


-7.8 


+7.5 


-6.9 


125.5 


19.15 


+7.2 


-7.8 


+7.5 


-6.9 


126.0 


19.01 


+7.2 


-7.8 


+7.5 


-6.9 


126.5 


18.86 


+7.2 


-7.8 


+7.5 


-6.9 


127.0 


18.71 


+7.1 


-7.8 


+7.5 


-6.9 


127.5 


18.57 


+7.1 


-7.8 


+7.5 


-6.9 


128.0 


18.43 


+7.1 


-7.8 


+7.5 


-6.9 


128.5 


18.29 


+7.1 


-7.8 


+7.5 


-6.9 


129.0 


18.15 


+7.1 


-7.8 


+7.5 


-6.9 


129.5 


18.01 


+7.1 


-7.8 


+7.5 


-6.9 


130.0 


17.88 


+7.1 


-7.7 


+7.5 


-6.9 


131.0 


17.62 


+7.1 


-7.7 


+7.5 


-7.0 


132.0 


17.36 


+7.0 


-7.7 


+7.5 


-7.0 


133.0 


17.11 


+7.0 


-7.7 


+7.4 


-7.0 


134.0 


16.86 


+7.0 


-7.7 


+7.4 


-7.0 


135.0 


16.62 


+7.0 


-7.7 


+7.4 


-7.0 


136.0 


16.38 


+6.9 


-7.6 


+7.4 


-7.0 


137.0 


16.14 


+6.9 


-7.6 


+7.4 


-7.0 


138.0 


15.92 


+6.9 


-7.6 


+7.4 


-6.9 


139.0 


15.69 


+6.9 


-7.6 


+7.4 


-6.9 


140.0 


15.48 


+6.9 


-7.6 


+7.4 


-6.9 


150.0 


13.53 


+6.7 


-7.4 


+7.4 


-7.0 


160.0 


11.85 


+6.5 


-7.3 


+7.5 


-7.1 


180.0 


8.810 


+6.2 


-7.0 


+7.4 


-7.5 


190.0 


7.815 


+6.1 


-6.9 


+7.4 


-7.5 


200.0 


7.082 


+6.0 


-6.8 


+7.4 


-7.7 


210.0 


6.496 


+6.0 


-6.7 


+7.4 


-7.8 


220.0 


6.005 


+5.9 


-6.6 


+7.3 


-7.6 



§ Notice that effectively the OFFP scheme matches the naive Breit Wigner implementation below 200 GeV. 
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m H (GeV) 


a (pb) 


scale(%) 


PDF+a s (%) 


230.0 


5.565 


+5.9 


-6.5 


+7.4 


-7.7 


240.0 


5.158 


+5.9 


-6.4 


+7.3 


-7.7 


250.0 


4.781 


+5.8 


-6.4 


+7.4 


-7.7 


260.0 


4.460 


+5.8 


-6.3 


+7.6 


-7.4 


270.0 


4.182 


+5.8 


-6.2 


+7.6 


-7.9 


280.0 


3.950 


+5.7 


-6.2 


+7.6 


-8.0 


290.0 


3.754 


+5.7 


-6.1 


+7.6 


-8.0 


300.0 


3.595 


+5.7 


-6.1 


+7.7 


-7.9 


310.0 


3.472 


+5.7 


-6.0 


+7.7 


-8.0 


320.0 


3.383 


+5.7 


-6.0 


+7.7 


-8.0 


330.0 


3.341 


+5.7 


-6.0 


+7.8 


-8.1 


340.0 


3.359 


+5.7 


-5.9 


+7.9 


-8.1 


350.0 


3.399 


+5.7 


-5.9 


+8.0 


-8.2 


360.0 


3.384 


+5.8 


-5.9 


+8.0 


-8.2 


370.0 


3.331 


+5.8 


-5.8 


+8.1 


-8.2 


380.0 


3.231 


+5.8 


-5.6 


+8.1 


-8.2 


390.0 


3.089 


+5.8 


-5.5 


+8.2 


-8.2 


400.0 


2.921 


+5.8 


-5.4 


+8.2 


-8.2 


420.0 


2.550 


+5.8 


-5.3 


+8.3 


-8.3 


440.0 


2.179 


+5.8 


-5.3 


+8.5 


-8.4 


450.0 


2.002 


+5.8 


-5.2 


+8.6 


-8.4 


460.0 


1.836 


+5.8 


-5.2 


+8.7 


-8.4 


480.0 


1.537 


+5.8 


-5.2 


+8.9 


-8.5 


500.0 


1.283 


+5.8 


-5.1 


+9.1 


-8.5 


520.0 


1.069 


+5.8 


-5.1 


+9.2 


-8.6 


540.0 


0.8911 


+5.8 


-5.1 


+9.4 


-8.6 


550.0 


0.8141 


+5.8 


-5.1 


+9.4 


-8.7 


560.0 


0.7442 


+5.9 


-5.1 


+9.4 


-8.7 


580.0 


0.6230 


+5.9 


-5.1 


+9.5 


-8.7 


600.0 


0.5231 


+5.9 


-5.0 


+9.5 


-8.8 


620.0 


0.4403 


+5.9 


-5.0 


+9.6 


-8.9 


640.0 


0.3719 


+5.9 


-5.0 


+9.7 


-9.0 


650.0 


0.3424 


+5.9 


-5.0 


+9.7 


-9.0 


660.0 


0.3153 


+5.9 


-5.1 


+9.8 


-9.1 


680.0 


0.2680 


+6.0 


-5.1 


+9.9 


-9.2 


700.0 


0.2289 


+6.0 


-5.1 


+10.1 


-9.3 


720.0 


0.1962 


+6.0 


-5.1 


+10.2 


-9.5 


740.0 


0.1687 


+6.1 


-5.1 


+10.4 


-9.6 


750.0 


0.1566 


+6.1 


-5.1 


+10.4 


-9.7 


760.0 


0.1455 


+6.1 


-5.2 


+10.5 


-9.7 


780.0 


0.1260 


+6.1 


-5.2 


+10.5 


-9.8 


800.0 


0.1095 


+6.1 


-5.2 


+10.6 


-9.8 



Table 1: Cross sections at the LHC (y/s = 8 TeV) and corresponding scale and PDF+ag uncer- 
tainties computed according to the PDF4LHC recommendation. 

The results we are going to present are obtained by using the MSTW2008 NNLO parton 
distribution functions (PDFs) [31], setting the reference values for the factorization and renor- 
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malization scales to the Higgs boson virtuality M 0. The on-shell width of the Higgs boson is 
evaluated with the program HDECAY [32J. We set the top-quark mass to m t = 172.5 GeV and 
we choose m b = 4.75 GeV and m c = 1.40 GeV consistently with the MSTW2008 set. 

Our predictions and the corresponding uncertainties, computed as discussed below, are pre- 
sented in Table 1. We stress that the inclusion of finite- width effects results in an increase of 
the cross section with respect to the ZWA of about 0(10%) at tjih = 800 GeV. It is well known 
that as mfl increases and finite-width effects become important, the signal cross section becomes 
itself ill defined, and for each decay channel of the Higgs boson, only the full signal + background 
computation in that channel strictly makes sensed. The use of a naive Breit Wigner, which would 
correspond to replace 7# with the on-shell width r(m#) in Eq. 01]), would give a smaller cross 
section with respect to the result in the complex-pole scheme, the difference ranging from —3.5% 
for m H = 300 GeV to -18% at m H = 600 GeV, to -27% at m H = 800 GeV. 

We now review the various sources of uncertainty affecting the cross sections presented in Table 
1. The uncertainty has two main origins: the one coming from the partonic cross sections, and 
the one arising from our limited knowledge of the PDFs. 

Uncalculated higher-order QCD radiative corrections are the most important source of uncer- 
tainty on the partonic cross section. A method, which is customarily used in perturbative QCD 
calculations, to estimate their size is to vary the renormalization and factorization scales around 
the hard scale M EE In general, this procedure can only give a lower limit on the true uncertainty. 
Here we quantify the uncertainty as in Refs. [21], [20]: we vary independently fip and [1r in the 
range 0.5M < Hf,H>r < 2M, with the constraint 0.5 < Hf/^r < 2. 

The scale uncertainty ranges from +7 — 8% (m# = 125 GeV) to about +6 — 5% (m# = 800 
GeV). The results are consistent with those of Ref. [21] [20]; in particular, we note that the 
effect of scale variations in our resummed calculation is considerably reduced with respect to the 
corresponding fixed-order NNLO result. 

Another source of perturbative uncertainty on the partonic cross sections comes from the 
implementation of the EW corrections. Our results are obtained in the complete factorization 
scheme, a scenario supported by the effective field theory computation of mixed QCD-EW cor- 
rections presented in Ref. [35]. The partial factorization scheme, in which EW corrections are 
applied only to the LO cross section, would lead to a change of our results varying from about 
—3% to +2% in the range of Higgs boson masses we consider. 

A potentially important source of perturbative uncertainty arises from the use of the large-m f 
approximation in the computation of the partonic cross section beyond NLL+NLO. The accuracy 
of the large-m t approximation at NNLO has been studied by computing subleading terms in the 
large-m 4 limit [SB], concluding that it works remarkably well, to better than 1% for mjj < 300 GeV. 
For heavier Higgs bosons, we expect the uncertainty due to the laige-mt approximation to play an 
increasing role. Nonetheless, it is well known that the comparison of the exact NLO calculation 
with the result obtained in the large mt limit but rescaled with the exact (m t and ra& dependent) 
LO cross section shows agreement within the 10% level even at high Higgs boson masses |37j . 

^The numerical integration over the virtuality M in Eq. (f4]) is performed between 50 to 1800 GeV. 
"A method to estimate the uncertainty from interference effects in the ZZ channel is proposed in Ref. 33 . 
**An attempt to go beyond this standard approach is made in Ref. [51] , 
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Since our calculation exactly includes the heavy-quark mass dependence up to NLL+NLO, the 
uncertainty due to the large-m t approximation should be well below the 10% effect in the mass 
range we consider. 

The other important source of uncertainty in the cross section is the one coming from PDFs. 
Our understanding of PDFs has improved considerably in the last years, and we have now various 
PDF sets at NNLO accuracy: MSTW2008 [3], NNPDF21 [38], JR09 [39] and ABK11 [40]. In 
order to produce the central values for the cross section, we rely on the MSTW2008 NNLO PDFs 
|31j . The PDF and as uncertainties are estimated using the corresponding 68% C.L. band from 
the MSTW sets, normalized according to the PDF4LHC prescription [41J. We note that the 
ensuing uncertainty is rather close to the one obtained by using the 90% C.L. set of MSTW. The 
uncertainty ranges from +8 — 7% {m H = 125 GeV) to +11 — 10% (m# = 800 GeV). By using the 
NNPDF21 NNLO default set [M] with a s (m z ) = 0.119 we find differences that range from +5% 
(m H = 125 GeV) to +2% (m H = 800 GeV) with respect to our central MSTW2008 result. When 
using the NNPDF21 NNLO set with as(mz) = 0.114 the difference ranges from —3% (m# = 125 
GeV) to -11% (m H = 800 GeV). 

The JR09 and ABM11 partons give larger differences with respect to our central MSTW 
prediction. For m H = 125 GeV the JR09 (ABM11) result is lower than MSTW08 by about 10% 
(7%). At larger Higgs masses the agreement of JR09 with MSTW08 improves, but the difference 
with ABM11 increases, being about 14% at m# = 300 GeV and further increasing at higher Higgs 
masses. We remind the reader that the ABM11 NNLO fit does not include Tevatron jet data and 
that the ensuing QCD coupling is as{rriz) = 0.1134, significantly smaller than the world average. 

Recently also the CTEQ collaboration has released an NNLO PDF set |44j . We find that the 
use of the CT10 NNLO central set, corresponding to asijnz) = 0.118 leads to results that agree 
at the 1% level or better with those obtained with MSTW2008. By using the sets corresponding 
to as = 0.116 and as = 0.120 the differences range from ±4% at m H = 125 GeV to +5 — 6% 
at = 800 GeV, thus well within our PDF4LHC uncertainty band. We find this agreement 
reassuring and we conclude that our central predictions, endowed with the PDF4LHC uncertainty, 
should be sufficiently conservative. Nonetheless, we believe that the large differences obtained with 
JR09 and ABM11 PDFs definitely deserve further investigations. 

The results of this paper can be compared to those presented in Ref. [19], where the impact of 
finite-width effects is studied as well. Besides the NLO QCD corrections with the exact dependence 
on the top- and bottom-quark masses, the NNLO corrections in the large- m t limit, and two loop 
EW effects [27], the calculation of Ref. [TU] includes mixed QCD-EW corrections evaluated in 
an effective field theory approach |35j and the independent evaluation of EW effects from real 
radiation [421 S3] whose effect is, however, at the 1% level or smaller, and that we neglect here. For 
a light Higgs boson, the main difference with our computation arises from the evaluation of higher- 
order QCD corrections. Following what previously done in Ref. [35], in Ref. [19] these corrections 
are computed up to NNLO but choosing \i F = fi R = m H /2, as an attempt to reproduce effects 
beyond NNLO, that, in our calculation, are instead estimated through soft-gluon resummation. 
For itih = 125 GeV the result of Ref. [19] with the corresponding scale uncertainty is a = 
20.69^g 3% pb, 7% higher with respect to our a = 19.31+2/g% pb, but still well within the uncertainty 
bands. Larger differences are observed in the high mass region, due to the different implementation 
of finite-width effects. In Ref. [19] a Breit-Wigner with running width is used as the default 
implementation of the line-shape. At m# = 400 GeV, the result of Ref. [T5] is about 16% smaller 
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than ours. 



In this paper we have presented updated predictions for the cross section for Higgs boson pro- 
duction at the LHC with y^s = 8 TeV, and discussed the corresponding uncertainties. The results 
are based on the most advanced theoretical information available at present for this observable, 
including soft-gluon resummation up to NNLL accuracy, two-loop EW corrections, exact treat- 
ment of heavy quark mass effects up to NLL+NLO accuracy and finite-width effects evaluated in 
the complex-pole scheme. We look forward to a comparison of our results with LHC data. 

We would like to thank Giampiero Passarino and Carlo Oleari for useful discussions. This 
work was supported in part by UBACYT, CONICET, ANPCyT and the Research Executive 
Agency (REA) of the European Union under the Grant Agreement number PITN-GA-20 10-264564 
(LHCPhenoNet). 
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